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ABSTRACT 


Kaib and Cowan (2015) recently used terrestrial planet formation simulations to conclude that the moon-forming 
irnpactor ( 'Lheia]~aad only a ^5% or less chance of having the same oxygen isotope composition as Earth, while 


Mastrobuono-Battisti et al. (2015) used seemingly similar simulations and methods to arrive at a highe r value of 
20% or mor e. Here we deriv e tne results of both papers from a single set of simulations. Compared to Kaib and 


Cowan| ( |2015[ ), the analysis of Mastrobuono-Battisti et al. (2015) systematically yields more massive Theia analogs 
and imposes flatter gradients across the original protoplanetary disk. Both of these effects diminish isotopic 

differences between Earth and Theia an alogs. While it is notori ously difficult to produce systems resembling our 


actual terrestrial planets, the analysis of 
analogs at orbital locations near the real p 

Kaib and Cowan 

(2015 

) more often selects and analyzes Earth and Mars 

anets. 

Given this, we conclude that the greater isotopic ditterences between 

Earth and Theia found in Kaib and Cowan 

(2015 

) better reflect the predictions of terrestrial planet formation models. 


Theia analogs consistent with the canonical giant impact hypothesis, this fraction still remains in the 5-8% range. 


1. BACKGROUND 

The canonical giant impact formation scenario for 
Earth’s moon predicts that a Mars-mass body, Theia, 
struck the proto-Earth in a glancing blow which threw 
material into orbit arou nd Earth and eventual l y co- 


alesced into the Moon (jHartmann and Davis 


1975 


Cameron and Ward 1976). Llydrodynamical modeling o: 


this collision has since shown that the Moon should con¬ 
sist of ^60-90% of material from Theia, with the rest 
comin g from the proto-Earth ( |Canup|[2004[ [Reufer et al. 
2012). Yet the Moon’s oxygen isotope composition is vir¬ 


tually identical to Earth (A^^O < 0.015%o) even though 
the isotopic compositions of nearly all meteorites differ 


significantly from the Earth ([Wiechert et al. 2001 Her- 


wartz et al. 2014|). In particular. Mars has a A^'O of 
0.32%o. Because simulations of terrestrial planet forma¬ 
tion predict that each large rocky body has a wide and 
stochastic feeding zone, it has been suggested that ter¬ 
restrial planet formation is unlikely to produce an Earth 
and Theia that are isotopically s imilar if other bodies (i.e. 
Mars ) are isotopically distinct (Pahlevan and Stevenson 
2007|. 


TCTrestrial planet formation simulations are an excel¬ 
lent tool to verify whether we should expect a large com¬ 
positional difference between Earth and its major im- 
pactors. These simulations model the assembly of sys¬ 
tems of terrestrial planets via the accretion of ^100-1000 


lunar-to-Mars size objects in orbit about the Sun (Ray¬ 
mond et al. 2014). Two recent independent works used 
suites of these simulations to statistically quantify com¬ 
positional differences between Earth analog p lanets and 
their last major impactors (Theia analogs) ( Kaib and 
Cowan 2015 Mastrobuono-Battisti et al. 2015 T While 
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very similar simulations were used, they arrived at signif¬ 
icantly different probabilities that Theia’s isotopic com- 
position was simi l ar to that of the Earth (and Moon). 


Kaib and Cowan (2015) (hereafter referred to as KC15) 
concluded that there was a <5% chance that Theia’s 
feeding zone was similar enough to the Earth’s to yield 
the M oon’s observed isotopic composition . On the other 
hand, Mastrobuono-Battisti et al. (2015) (hereafter re¬ 
ferred TO’as^MPHTFjTounTa'Ti^ablllfy of >20%. 

In this brief follow-up work, we report the cause of 
this discrepancy. The main differences between the as¬ 
sumptions used in KC15 and MPR15 are: (1) the criteria 
employed to select Earth and Mars analog pairs used to 
calibrate the initial A^^O distribution in each simulation, 
and (2) the criteria employed to select Earth and Theia 
analog pairs used to determine the oxygen isotopic dif¬ 
ference between Earth and Theia in each simulation. In 
KC15, Earth analogs are required to have masses over 
0.5 Mq and orbits between 0.8 and 1.2 AU, while Mars 
analogs are assumed to be the next planet outward with a 
mass greater than 0.05 M0. Once the Earth-Mars ana¬ 
log pair is selected, the accretion history of the inner 
planet (the Earth analog) is searched for an irnpactor 
with a mass larger than 0.1 M0. If such impactors exist, 
the last of these is designated a Theia analog, and this 
body’s feeding zone is compared against the Earth ana¬ 
log’s. The isotopic composition of the Earth and Theia 
analogs is then calculated by imposing a A^^O gradi¬ 
ent on the original planetesimals that yields a A^'O of 
0.32%o between the Earth analog and the Mars analog. 

Meanwhile, MPR15 place more emphasis on planet or¬ 
dering than orbital elements or planet mass to select 
Earth-Mars analogs. In systems that form 4 or more 
planets, the 3rd and 4th planets are assumed to represent 
Earth and Mars, respectively, while in 3-planet systems, 
Earth and Mars are taken to be the 2nd and 3rd planets 
out from the Sun. These Earth-Mars pairs are selected 
regardless of their masses or orbits. These planet pairs 
are then used to set the primordial A^^O gradient just 
as in KC15. Next, any planet in the system is searched 
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Fig. 1.— The distribution of |A^^O| for Theia analogs when 
employing the Earth, Mars, and Theia analog selection criteria of 
MPR15 (dashed blue) and when employing the analog selection 
of KC15 (red solid). The criteria are applied to the same sets of 
simulations (CJS and EJS from KC15). The shaded region marks 
Theia analogs with |A^^O| values that are less than or equal to 
that of the Moon (0.015 %o). 

for impactors composed of at least 50 particles, and these 
are taken to be Theia analogs. These impactors need not 
strike the Earth analog used to set the isotopic gradient 
and typically strike another planet in the system instead. 
Finally, the isotopic composition of the Theia analog and 
its target planet are compared by employing the isotopic 
gradient imposed by the Earth-Mars analog pair. 


2. ANALYSIS AND RESULTS 

To compare the KC15 and MPR15 analyses, we use 
the 100 CJS and EJS simulations from KC15 since these 
are virtually equivalent to those employed in MPR15. 
(It should be noted that the EJS initial conditions are 
effectively equiv a lent t o the EEJS initial conditions in 
Raymond et al. (2009).) KC15 contains a detailed de- 
scription of the simulations. In KC15, the analysis of 
CJS and EJS simulations were always performed sepa¬ 
rately to control for any effect of different giant planet 
configurations, while MPR15 combined numerous sets of 
simulations that used a variety of giant planet configu¬ 
rations. Because there was no distinction made between 
initial conditions sets in MPR15, we combine the CJS 
and EJS simulations for this follow-up work. 

We now apply the analysis described in MPRI5 to the 
combined CJS and EJS simulations to derive the iso¬ 
topic composition of Theia analogs. The results of this 
analysis are shown in Figure I. Here we plot the dis¬ 
tribution of |A^^O| values for Theia analogs relative to 
the impacted planet in these simulations. Our modified 
analysis predicts that 15% (10 out of 65) of our Theia 
analogs have |A^'^0| < 0.015%o, very near the results of 
MPR15. Recall that the KC15 results derived from the 
same simulation sets predict a probability of ~2% (1 out 
of 53) that Theia has an Earth-like composition. Thus, 
we recover the results of both KC15 and MPR15 from 
the same set of simulations. 

Several different factors cause the differing isotopic dis¬ 
tributions for Theia in Figure 1. First, MPR15 require 
impactors (Theia analogs) to be comprised of at least 50 
original particles, whereas KC15 only impose a minimum 
Theia mass of 0.1 M® and their median Theia analog is 
comprised of 33 particles. While the 50-particle limit 


minimizes mass resolution effects, it also causes the im- 
pactor masses to be larger in MPR15 than in KC15. In 
the giant impact scenario, the proto-Earth still accretes 
nearly all of the impactor, and a more massive impactor 
pulls the final isotopic composition of the Earth closer 
to that of Theia. Indeed, when MPR15 relax their 50- 
particle requirement in their Extended Data Table 2, a 
clear trend toward more isotopically dissimilar impactors 
is seen. (In the hit-and-run M oon formation scenario 
proposed by Reufer et al. (2012), the Earth does not ac¬ 
tually accrete most of the impactor, and the impactor’s 
influence on the Earth’s final composition would be di¬ 
minished. However, neither MPR15 or KC15 considered 
this complication.) 

When we perform the MPR15 analysis on our simu¬ 
lations, the median mass ratio of Theia to the proto- 
Earth is 0.35. Meanwhile, the KC15 analysis yields a 
median mass ratio of 0.25, and only 17% of KC15 Theia 
analogs exceed mass ratios of 0.35. Most impacts involv¬ 
ing such large mass ratios are typically found for planets 
well inside 1 AU, nearer the inner disk edge where par¬ 
ticle density and accretion rates are higher. Indeed, the 
planetary targets of the moon-forming collisions selected 
via the MPR15 criteria have a median semimajor axis 
of 0.75 AU. This places many target planets near the 
hard inner disk edge imposed by initial conditions, and 
these planets may also have inherently different feeding 
zones compared to the more distant target planets (Earth 
analogs) required in KC15. 

In addition, because the isotopic difference between the 
Earth and Mars analog pair is always fixed at 0.32%o, the 
A^^O gradient imposed on the initial disk of planetesi- 
mals is completely dependent on the separation between 
the mean formation distance of the Mars analog and the 
mean formation distance of the Earth analog. A larger 
separation will produce a flatter gradient, which will in 
turn yield a smaller difference between the Earth’s and 
Theia’s isotopic composition (they are identical if the 
gradient is completely flat). In many instances, the 4th 
planet of simulated systems (the Mars analog under the 
MPRI5 criteria) orbits well beyond 2 AU. This often 
leads to orbital separations of over ~I AU between the 
Earth and Mars analogs used to set the A^^O gradient 
in MPR15. 

If these widely separated planets also have widely sep¬ 
arated mean formation distances, we might expect that 
some of the success cases of MPR15 come from systems 
with unusually flat A^^O gradients. In fact, MPR15 
explored the effects of employing steeper A^^O gradi¬ 
ents in their Extended Data Table 2 and found a signif¬ 
icant decrease in the fraction of isotopically Earth-like 
Theia analogs if the disk A^^O gradients were increased 
by 25%. In addition, we find a correlation coefficient of 
0.72 when we compare the mean formation separations 
with the final orbital separations for the MPRI5 Earth- 
Mars pairs. Given this, we may expect that the widely 
separated Earth-Mars pairs of MPRI5 will in fact have 
widely separated mean formation distances, resulting in 
flat A^’^0 gradients. 

The separations of the mean formation distances for 
Earth-Mars pairs selected with the KC15 and MPR15 
criteria are shown in Figure 2. As expected, we see that 
the Earth-Mars pairs selected by the MPRI5 prescription 
have a long tail to very high separations not seen among 
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Fig. 2. — The distribution of separations between the mean for¬ 
mation radii of Mars analogs and Earth analogs as selected by the 
criteria of MPR15 {dashed blue) and by the criteria of KC15 {red 
solid). The distribution of Earth-Mars separations is also shown for 
simulations that yielded Theia analogs with |A^^O| values equal to 
or less than that of the Moon when employing the MPR15 criteria 
{dotted black). 

the Earth-Mars pairs selected with the KC15 criteria. As 
mentioned previously, performing the MPR15 analysis on 
our simulations yields 10 instances where Theia’s |A^^O| 
is less than or equal to the Moon’s. As Figure 2 shows, 6 
of these 10 instances have Earth-Mars pairs that lie in the 
tail of very high Earth-Mars separations, which require 
flat A^’^0 gradients. This suggests that the different dis¬ 
tributions of Earth-Mars analog formation distances is 
a primary reason for the increased likelihood that Theia 
has an Earth-like isotopic composition in MPR15. 

The large separations between the mean formation dis¬ 
tances of Earth analogs and Mars analogs seem to only 
occur in systems where at least one of the two analog 
planets has a very small mass. In Figure 3 we plot the for¬ 
mation distance separation against the minimum planet 
mass of each Earth-Mars analog pair. Here we see that 
the separation between Earth’s and Mars’ formation dis¬ 
tance is only above 1 AU if at least one of the analog 
planets has a mass well below 0.1 M 0 . These very low 
mass planets are comprised of one to a few particles, and 
it has been shown that individual particles can migrate 
very large d istances during terrestrial planet formation 


simulations i 

Chambers|2001 |Raymond et al.|2004 Mor- 

ishima et al. 

2U1UI. This causes the large formation dis- 


tance separations and the flat A^^O gradients that are 
consequently imposed. If we only consider cases where 
both Earth and Mars analogs have masses above 0.06 
Mq (approximately the smallest planetary mass seen in 
our Solar System) then the fraction of MPR15 Theia 
analogs with |A^’^0| < 0.015%o drops to just 3% (1 out 
of 31 cases), nearly the same as the KC15 results. 

It should be noted, however, that Mars analogs should 
in fact be substantially smaller than the Earth analogs, 
and no such requirement was made in KC15. If we use 
the KC15 Earth-Mars pairs and now require that the 
Earth analog is at least twice as massive as the Mars 
analog, our usable sample of simulations drops from 53 
to 18. Of these 18 systems, we still only have one Theia 
analog with | A^^O| < 0.015%o (a 5.5% success rate), and 
the median value of |A^^O| actually increases slightly 
(0.16%o vs. 0.18%o). Thus, the results of KC15 do not 
seem to be due to having selected Mars analogs that are 



(Mg) 

Fig. 3. — The separation between the mean formation distances 
of Mars analogs and Earth analogs as selected by the criteria of 
MPR15 {dashed blue) are plotted against the minimum planet mass 
of each Earth-Mars pair. The shaded region highlights systems 
where the mass of Earth and/or Mars is below 0.06 M0. 

unrealistically massive. 

2 . 1 . Simulation Uncertainties 

The stochastic nature of planet formation is well cap¬ 
tured by numerical experiments: if a simulation is re¬ 
peated with virtually identical initial conditions, plan¬ 
ets form at different locations and have different feed¬ 
ing zones. This introduces an uncertainty in the results 
of any individual simulation, and is precisely the reason 
that KCI5 and MPRI5 ran large ensembles of simula¬ 
tions to quantify the likelihood of the Moon-forming im¬ 
pact. 

The “granular” nature of N-body simulations also in¬ 
troduces an uncertainty in the planetary composition. 
Our simulations contain 1000 planetesimals to represent 
small bodies, but the real number may have been orders 
of magnitude higher during Solar System formation. The 
feeding zone of each simulated protoplanet is likely con¬ 
structed from a much smaller number of planetesimals 
than in the actual solar system. 

Granularity is a consequence of computational limi¬ 
tations and can artificially skew our results. Even if the 
true feeding zones of a Theia analog and an Earth analog 
were identical, our simulations could result in a non-zero 
|A^'^ 0 | due to the small number of simulated planetes¬ 
imals. Alternatively, the coarse mass resolution of our 
simulations could produce two planets with nearly iden¬ 
tical compositions where a higher-resolution simulation 
would have shown that they have significantly different 
feeding zones. 

Disentangling granularity error from natural stochas- 
ticity is difficult. MPR15 estimate the uncertainty in 
each simulated major body’s mean formation distance 
using the standard error of the mean (SEM). The SEM 
values for Earth and Theia analogs are then added in 
quadrature and translated to a A^"^© uncertainty. When 
simulated | A^^O| values are below O.OI5%o-l- the SEM un¬ 
certainty, MPR15 deem that the simulation successfully 
explained the Moon’s isotopic composition. This only 
makes sense, however, if one presumes that the feeding 
zones of Earth and Theia are inherently identical and 
that simulation granularity artificially increased compo¬ 
sitional differences. 

We instead estimate the granularity error of our sim- 
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ulations by resampling each major body’s planetesimal 
distribution via bootstrap. If a simulated planet accreted 
N planetesimals from different semi-major axes, then we 
draw N times from that planetesimal distribution, with 
the possibility of repeated selections. This produces a 
new distribution of N planetesimals. By repeating this 
process 1000 times for each body, we generate a distri¬ 
bution of feeding zone locations for each simulated body. 
Bodies made of a few planetesimals from disparate semi¬ 
major axes will have higher variance than bodies built 
from many nearby planetesimals. The standard devi¬ 
ation of the feeding zone location for each body is an 
estimate of the granularity uncertainty, and is typically 
within a few percent of the SEM value. 

We only perform bootstrapping and SEM for distribu¬ 
tions of planetesimals, which are all equal-mass particles. 
The number of simulated embryos is comparable to ex¬ 
pected number in the real solar system, so our embryo 
population should not suffer from the granularity issues. 
It is not clear whether or not MPR15 make a distinc¬ 
tion between embryos and planetesimals in their SEM 
calculations. 

We use bootstrap resampling to predict—for each 
simulation—the probability that the feeding zones yield 
a Theia analog with |A^^O| < 0.015%o. For each in¬ 
dividual set of Earth-Mars-Theia feeding zones, we use 
bootstrap resampling to generate a distribution of 1000 
gradients and A^^O values for each individual sim¬ 
ulation. We then sum up the total probability that 
|A^’^0| < 0.015%c for our ensemble of simulated Theia 
analogs. Unlike the raw simulation results, these re¬ 
sults will account for the coarse mass resolution of our 
simulations. This analysis predicts a ~4.8% probability 
that Theia will have |A^^O| < 0.015%c. This number is 
2.5 X larger than the predictions of our raw simulation 
results, but still significantly smaller than the prediction 
of MPR15. 

Finally, it is worth noting that both SEM and boot¬ 
strap resampling implicitly assume that each accreted 
object represents an independent sample of a proto¬ 
planet’s “true” underlying feeding zone. During accre¬ 
tion, however, small bodies often merge with each other 
before being accreted onto even larger bodies. Since 
much of a protoplanet’s feeding zone is derived from con¬ 
glomerations of associated bodies, the initial locations of 
accreted planetesimals are likely correlated. It is there¬ 
fore possible that neither SEM nor bootstrapping prop¬ 
erly quantifies the uncertainty arising from using too few 
planetesimals. 


2.2. Earth’s Contribution to the Moon 

In KC15, it was assumed that lunar material is almost 
entirely derived from Theia instead of the Earth. In re¬ 
ality, Earth always contributes some fraction of material 
toward the Moon’s assembly. In the canonical giant im¬ 
pact hyp othesis, the E arth fraction is typically around 


15-25% (Canup 2004). In alternative scenarios, such 
as hit-and-run collisions, high velocity impacts, or very 


material can exceed 40% ( 

Reufer et al. 

2012 

Cuk and 

Stewart 

2012 

Canup 2012 

). Considering a significant 


analogs with larger |A^'0| values to explain the Moon’s 



Fig. 4.— The probability that Theia will have a value capa¬ 

ble of producing a moon with |A^^O| < 0.015%o is plotted against 
the assumed fraction of lunar material derived from the Earth. 
This probability is calculated from our bootstrapped A^^O {solid 
black) distribution, which accounts for granularity uncertainty, and 
from our raw simulation results {dotted black). 

isotopic composition since the | A^^O| will be scaled down 
by the compositional fraction that Theia provides (1 mi¬ 
nus the terrestrial fraction) when Theia material is in¬ 
corporated into the Moon. In Figure 4, we examine how 
the fraction of Theia analogs that successfully explain the 
lunar isotopic composition varies as we consider greater 
contributions from the Earth. Using our |A^^O| values 
generated from bootstrap resampling, we find that our 
Theia success rate moves from 4.8% to 8.4% as we con¬ 
sider Earth contributions up to 40% of the Moon’s mass. 
When looking at the raw simulation data instead of our 
resamplings, one finds an increase from ~2% with no ter¬ 
restrial contribution up to ^9.5% if the Earth supplies 
40% of the Moon’s material. 

Unlike KC15, MPR15 did consider scenarios in which 
the Moon had a significant terrestrial component and 
found that a terrestrial contribution of 40% could in¬ 
crease the success rate of Theia analogs up to 25% or 
55%, depending on whether their SEM-derived uncer¬ 
tainties were being considered. As shown in Figure 4C, 
KC15 finds success rates that are significantly lower than 
either of these values. We attribute this difference to 
the Earth, Mars, and Theia analog selection criteria em¬ 
ployed in KC15 vs MPR15. 


3. DISCUSSION AND CONCLUSIONS 

Although recent terrestrial planet formation models 
have shown more promise in producin g Martian analogs 
(Walsh et al.|2011 Izidoro et al.|2014 1, many of the sim- 
ulations employed in KC15 ancf 1 V1FR15 have difficulty 
reproducing the low mass of Mars (Raymond et al. 2009). 
Because of this, it is not obvious how to designate Earth 
and Mars analogs. However, we believe the KC15 crite¬ 
ria are more self-consistent for several reasons. First, the 
target of the moon-forming collision is always required to 
be the same planet used in calibrating the initial A^^O of 
simulation particles, whereas the target planets selected 
under the criteria of MPR15 often orbit inside Venus’ 
orbit. In addition, to minimize simulation resolution ef¬ 
fects, MPR15 require Theia analogs to be accreted from 
at least 50 particles, while the median Theia analog in 
KC15 is comprised of 33 particles. However, this 50- 
particle limit forces Theia analogs to always have larger 
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masses than that predicted in the canonical giant impact 
hypothesis. Finally, although the masses of Mars analogs 
are often too large in KC15, all of the simulations that 
are used in this work have a planet orbiting near 1.5 
AU, just as the Solar System does, and this is used to 
calibrate the initial distribution. Meanwhile, the 

Earth-Mars pairs selected with the MPR15 criteria can 
often orbit in the 2-3 AU region and have masses similar 
to or less than Mercury, something unseen in our Solar 
System. Flawed as these simulations may be at repli¬ 
cating the detailed structure of the inner Solar System, 
they do predict that the last major impactor on a mas¬ 
sive planet near 1 AU is quite unlikely to have the same 
isotopic composition as the planet itself, when there is 
also a planetary mass body near 1.5 AU with a Mars- 
like composition. 

MPR15 correctly point out that the feeding zones de¬ 
rived from terrestrial planet formation simulations are 
artificially influenced by the mass resolution of Simula- 
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tions. This increases the uncertainty in feeding zone lo- 
cations, and hence values. When we account for 

this uncertainty, the fraction of Theia analogs capable of 
explaining the Moon’s isotopic similarity to the Earth in¬ 
creases from 2% to 4.8%. Accounting for the possibility 
that the Moon could contain up to 40% terrestrial ma¬ 
terial further increases the successful fraction of Theia 
analogs to 8.4%. However, this is still well below the 
~50% success rates reported by MPR15. Given this, we 
conclude that standard models of the final assembly of 
terrestrial planets predict there is a ~5% chance that 
the canonical giant impact hypothesis can account for 
the isotopic composition of the Moon. 
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